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ABSTRACT 

We investigate equilibrium sequences of relativistic stars containing purely toroidal magnetic fields 
with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), EPS (Pandharipande et 
al.), Shen (Shen et al.), and LS (Lattimer & Swesty). We numerically construct thousands of equilib- 
rium configurations in order to study the effects of the realistic EOSs. Particularly we pay attention 
to the equilibrium sequences of constant baryon mass and/or constant magnetic flux, which model 
evolutions of an isolated neutron star, losing angular momentum via gravitational waves. Important 
properties obtained in this study are summarized as follows ; (1) Unlike the polytropic EOS, it is found 
for the realistic EOSs that the maximum masses do not monotonically increase with the field strength 
along the constant magnetic flux sequences. (2) The dependence of the mass-shedding angular velocity 
on the EOSs is determined from that of the non-magnetized case. The stars with Shen(FPS) EOS 
reach the mass-shedding limit at the smallest (largest) angular velocity, while the stars with SLy or 
Lattimer-Swesty EOSs take the moderate values. (3) For the supramassive sequences, the equilibrium 
configurations are found to be generally oblate for the realistic EOSs in sharp contrast to the poly- 
tropic stars. For FPS(LS) EOS, the parameter region which permits the prolately deformed stars is 
widest (narrowest). For SLy and Shen EOS, it is in medium. Furthermore, the angular velocities fiup, 
above which the stars start to spin up as they lose angular momentum, are found to depend sharply 
on the realistic EOSs. Our analysis indicates that the hierarchy of this spin up angular velocity is 
^^up.SLy > ^^up,FPS > ^^up,LS > ^^up.shon and this relation holds even if the sequences have strong 
magnetic fields. Our results suggest the relativistic stars containing purely toroidal magnetic fields 
will be a potential source of gravitational waves and the EOSs within such stars can be constrained 
by observing the angular velocity, the gravitational wave, and the signature of the spin up. 
Subject headings: magnetic fields, relativity, stars: neutron and magnetic fields 

1. INTRODUCTION 

Neutron stars observed in nature are magnetized with the typical magnetic field strength ~ 10^^-10^'^ G 
()Lvne fc Graham-SmithI [20051 ). For a special class of the neutron stars such as soft gamma-ray repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs), the field strength is often much larger than the canonical value as ~ 10 ^^ G, 
and these objects are collectively referred as magnetars (jLattimer &: Prakashll2007t IWoods fc Tho mpson"2004'). Al- 
though such stars are estimated to be only a subclass (~ 10%) of the canonical neutron stars ([Kouveliotou 1998), 
much attention has been drawn because they pose many astrophysically exciting but unresolved problems. 

Giant flaring activities observed in the SGRs have given us good o pportunities to study the coupling of the interior 
to the magnetospheric structures ([Thompson fc Duncan|[l995l |1996[ ). but we still kn o w litt le of the relationship be- 
tween the crustal fraction and the subsequent starquakes (see references in IWattj ()2006() : iGeppert fc Rheinhardtj 
(j2006i)). The origin of the large magnetic fi eld is also a big problem, whether descended from the main se- 
quence stars (iFerrario fc Wickramasingh j l2006f ) or generated at post-collapse in the rapidly rotating neutron star 
([Thompson fc Duncanlll993[ ). Assuming large magnetic fields before core-collapse, extensive m agnet ohydrodynamic 
(MHD) stellar collapse simulations have been carried out recentlv ([Kotake et al.l [2004t [Obc rgaulin ger et all 20061 : 
Shibata et al."2006 ULivne et al.l[2007t [Dessart et aD[2007[ : iTakiwaki et al.l[2004[ : [Sawai et al.[[2008t iKiuchi et al.[ [20081 : 
Takiwaki et al. 200^) towards the understanding of the formation mechanism of magnetars. Here it is worth mention- 
ing that the gravitational waves could be a useful tool to supply us with the information about magnetar interiors 
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(jBonazzola fc Gourgoulhoiilll996l : ICutleij|2002f ). While in a microscopic point of view, effects of ma gnetic fields larger 
than the so-called QED limit of i?QED = 4.4 x 10^^ G, on the EOSs (e . g-.lLattimer fc Prakashl (|2007| )) and the radiation 
processes have been also elaborately investigated (see lHarding fc Lai ( 20061 ) for a review). For the understanding of the 
formation and evolution of the magnetars. the unification of these macroscopic and microscopic studies is necessary, 
albeit not an easy task. 

In order to investigate those fascinating issues, the construction of the equilibrium configuration of magnet a rs ma y 
be one of the most fundamental problems. Starting from the pioneering study bv IChandrasekhar fc Ferinil ( 1953I) . 
extensive s tudies have bee n done and th is resea rch field is now experiencing "renaissance" (iBoucauet et al. 19951: 



loka fc Sasaki 2003. 2004; Kiuchi fc Kotak c 2008; iKonno et aLlll999l : iTomimura fc Eriguchill2005t lYoshida fc Eriguchil 



20061: lYoshida et al.ll2006t ). in which different levels of sophistication in the treatment of the magnetic field structure. 



equations of state (EOSs), and the general relativity, have been undertaken. Among them, a more sophistication is 
required for the studies employing the Newtonian gravity, because the general relativity(GR) should play an important 
role for the equilibrium configurations of compact objects like neutron stars. As mentioned below, typical densities 
for neutron stars interior exceed the nuclear density ~ lO^^g/cm"^. In such a high density region, the pressure P and 
rest mass density po becomes comparable, namely P ~ poc with c being the speed of light. In such a regime, the 
Newtonian gravity is too weak to gravitationally bind the neutron stars without the general relativity, that is GM/c^R 
becomes the orders of magnitude 10~^ with G, M and R being gravitational constant, mass and radius of star. 

It is noted here that the equilibr ium configurations of r elativistic stars without magnetic fields have been elaborately 
studied using the LORENE code (jBonazzola et al.iri998[ ). Unfortunately, however, the method of constructing a fully 
general relativistic star with arbitrarily magnetic structures is still not available . In most previou s stu dies, the weak 
magne tic fields and/or purely poloidal magnetic fields have been assumed. Boucquet et al.l (Il995l) and ICardall et aTl 
(|200lh have treated relativistic stellar models containin g purely poloidal ma gnetic fields. iKonno et al.l ( 19991 ) have 
analyzed similar models using a perturbative approach. Iloka fc Sasakil (j2004[ ) have investigated structures of mixed 
poloidal-toroidal magnetic fields around a spherical star by using a p erturbative technique . 

As shown in the MHD simulations of core-collapse supernovae (see iKotake et al.l ()2006[ ) for a review) , the toroidal 
magnetic fields can be efficiently amplified due to the winding of the initial seed poloidal fields as long as the core 
rotates differentially. After core bounce, as a result, the toroidal fields generically dominate over the poloidal ones even 
if there is no toroidal field initially. Eve n for the weakly magnet ized star prior to core-collapse, it should be mentioned 
that the magnetorotational i nstabilitv (iBalbus fc Hawlevlll99l[ l co uld be anoth er effective mechanism for generating 
large toroidal magnetic fields (jAkivama et al.l[2"003f ). As pointed bv I Spruit (|2007f ). such large magnetic fields generated 
by those processes may be confined inside the neutron star for thousands of years. Moreover, recent calculations of 
stellar evolution suggest that toroidal component of magnetic field domin ates over poloidal one with the orders of 
magnitude 10~^ in the newly born neutron star after supernova explosion (jHeger et al.ll2005l ). These outcomes seem 
to indicate that some ne utron stars could have toroidal fields much higher than the poloidal ones, which motivated 
iKiuchi fc Yoshidal ()2008D to study the equilibrium configurations with purely toroidal fields. In the study, the master 
equations were derived for the first time for obtaining the relativistic rotating st ars with the t oroidal fields. On the 
other hand, it is well known that pure toroidal magnetic field in star is unstable (|Tavleilll973f ). iBraithwaite fc SpruitI 
(|2004f ) have investigated the evolution of magnetic field in stellar interior and found "twisted torus" field configuration 
is stab le, in which the poloidal and toroidal magnetic field have comparable field strength. Recently, iKiuchi et al.l 
(|2008f l have confirmed the Tayler instability in the neutron stars with the general relativistic magnetohydrodynamics 
simulation (GRMHD) under axisymmetric condition. Important finding of this study is the toroidal magnetic fields 
settle down to the new equilibrium states with the circular motion in meridian plane. This suggests that the toroidal 
magnetic field may remain in neutron stars even if the poloid al magnetic field is much weaker than the toroidal one. 

However there may remain a room yet to be sophisticated in IKiuchi fc Yoshidal (|2008[ ) that the polytropic EOSs have 
been used for simpli city. In general, the c entral density of neutron stars is considered to be higher than the nuclear 
saturation density ( Bavm fc Pethiclall979D . Since we still do not have a n definite answer about the EOS in such a 
higher regime, many kinds of the nuclear EOSs have been proposed (e.g.. iLattimer fc Prakaslil (120071)') depending o n 
the descriptions of the nuclear forces and the possible appearance of the exotic physics (e.g., iGlendennind (|200lh ). 
While the stiflFness of the polytropic EOS is kept constant globally inside the star, it should depend on the density in the 
realistic EOSs. Since the equilibrium configurations are achieved by the subtle and local balance of the gravitational 
force, centrifugal force, the Lorentz force, and the pressure gradient, it is a nontrivial problem how the equilibrium 
configurations change for the realistic EOSs. 

These situations motivate us to study the effects of the realistic EOSs on the general relativistic stellar equilib- 
rium configurations with purely toroidal m a gnetic fields. For the purpose, we in corporate the realistic EO Ss to 
the numerical scheme bv IKiuchi fc Yoshidal (I2008D. Four ki nds of EOSs, of SLy (iDouchin fc Haensell l200lh FPS 
(jPandharipande fc Ravenhalllll989D . Shen (jShen et all (199831) . and Lattimer-Swestv(LS) (jLattimer fc Swestvl [l99ll) 
are adopted, which are often employed in the recent MHD studies relevant for magnetars. By so doing, we con- 
struct thousands of equ ilibrium stars for a wide range of parameters with different kinds of EOSs. As shown by 
iTrehan fc Uberoil (|1972D . the toroidal field tends to distort a stellar shape prolately because the toroidal field lines 
behave like a rubber belt, pulling in the matter around the magnetic axis. Such a prolate-shaped n eutron star i s 
pointed out to be optimal gravitational- waves emitters, depending on the degree of the prolateness ()Cutleill200l l. 
Thus, we pay attention to how the prolateness is affected by the realistic EOSs. For making such configurations, 
maximum magnetic field strength deep inside the magnetars c an become as high as 10^^ G, whic h could be generated 
due to a — dynamos in the rapidly rotating neutron stars ((Thompson fc Duncanl[l993l Il996l ). It should be noted 
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that the possibility of such ultra-magnetic fields has not been rejected so far, because what we can learn from the 
magnetar's observations by their periods and spin-down rates is only their surface fields 10^^ G) . We construct 
equilibrium sequences, along which the total baryon rest mass and/or the magnetic flux with respect to the meridional 
cross-section keep constant. It has been s uggested that the magnetic field will decay by the process of the Ohmic decay, 
ambipolar diffusion, and Hall drift (e.g.. iGoldreich fc Reisenegged ()1992l )). This implies the magnetic flux might not 
be well conserved during quasi-steady evolution of isolated neutron stars. However, taking account into these effect 
to the neutron star's evolution requires extensive efforts. Hence, as a first approximation, we explore the evolutionary 
sequences by use of the equilibrium sequences of constant rest mass and magnetic flux. 

This paper is organized as follows. The master equations for obtaining equilibrium conflgurations of rotating stars 
with purely toroidal magnetic fields are briefly reviewed in Sec. [51 The numerical scheme for computing the stellar 
models and EOSs we employ in this study are briefly described in Sec. [3l Sec. |4] is devoted to showing numerical 
results. Summary and discussion follow in Sec. [3 In this paper, we use geometrical units with G = c = 1. 

2. SUMMARY OF BASIC EQUATIONS 

Master e quations for the rotating relativistic stars containing purely toroidal magnetic flelds have been derived in 
iKiuchi (SpY oshida (2008). Hence, we only give a brief summary for later convenience. 

Assumptions to obtain the equilibrium models are summarized as follows ; (1) Equilibrium models are stationary 
and axisymmetric. (2) The matter source is approximated by a perfect fluid with infinite conductivity. (3) There is 
no meridional flow of the matter. (4) The equation of state for the matter is barotropic. Although we employ realistic 
EOSs in this paper, this barotropic condition can be maintained as explained in subsection 13. 21 (5) The magnetic axis 
and rotation axis are aligned. 

Becaus e the circularity condi tion ( see, e.g. [W ald (1984')) holds under these assumptions, the metric can be written, 
following iKomatsu et al.l (|1989( 1 and lCook et al.. (1992 ). in the form, 

ds^ = -e^+Pdt^ + e^^idr^ + T^d0^) + e'^^V sin^ 9{dip - udtf , (1) 

where the metric potentials, 7, p, a, and uj, are functions of r and 6 only. We see that the non-zero component of 
Faraday tensor F^i, in this coordinate is Fi2. In case of purely toroidal magnetic held, the Grad-Shafranov equation is 
known to be trivial because it is given only in terms of the (p component of the vector potential (see Kiuchi & Yoshida 
(2008) for details). Then the toroidal magnetic field can be determined by the integrability condition of the equation 
of motion (see Eq.(3)). Integrability of the equation of motion of the matter requires, 

e'^-2a sin6IFi2 = K{u); u = pf^he^^r^ sin^ 6*, (2) 

where K is an arbitrary function of pohe^^r^ sin^ 9. The variables po and h represent the rest mass density and 
relativistic specific enthalpy, respectively. Integrating the equation of motion of the matter, we arrive at the equation 
of hydrostatic equilibrium, 

where v = {VL — iS)r sm9e~P with Q. being the angular velocity of the matter and C is an integration constant. 

Here, we have further assumed the rigid rotation. To compute specific models of the magnetized sta rs, we need specify 
the fu nction forms of K, which determines the distribution of the magnetic fields. According to iKiuchi fc Yoshidal 
(|2008l ). we take the following simple form, 

K{u) = bu\ (4) 

where b and k are constants. Regularity of toroidal magnetic fields o n the magnetic axis requi res that fc > 1. If fc > 1, 
the magnetic fields automatically vanish at the surface of the star. IKiuchi fc Yoshidal ()2008f ) have investigated how 
the choice of fc affects the mag netic field d i stribu tion. In this study, we consider the fc = 1 case because in the general 
relativistic MHD simulation, K iuchi et al.l (l2008l) have found that magnetic distribution with fc 7^ 1 is unstable against 
axisymmetric perturbations. Note that lTavleil ()1973l ) has suggested that the magnetic field configuration of fc = 1 is 
unstable against non-axisymmetric perturbations. 

3. NUMERICAL METHOD AND EQUATION OF STATE 

3.1. Numerical method 

To solve the master equations numerically, we employ the Kiuchi- Yoshida scheme (jKiuchi fc Yoshidal l2008l ). which 
does not care about the function form of EOS. Hence, it is straight forward to update our numerical code for incorpo- 
rating the realistic EOSs discussed below. 

After obtaining solutions, it is useful to compute global physical quantities characterizing the equilibrium configu- 
rations to clearly understand the properties of the sequences of the equilibrium models. In this paper, we compute 
the following quantities: the gravitational mass M, the baryon rest mass Mq, the total angular momentum J, the 
total rotational energy T, the total magnetic energy H, the m agnetic flux ^, the gravi tational energy W and the 
mean deformation rate e, whose definitions are explicitly given in IKiuchi fc Yoshidal ()2008f ). More explicitly, the mean 
deformation rate e is defined as 
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TT J er"* sin 0(1 + cos^ 6)drd9 and Izz = 27r J er^ sin'^ OdrdO with e being the energy density of the matter, 
(g) Circumferential radius i?cir is defined as i?cir = c'^^^'^^^'^re with re being the coordinate radius at the stellar 
equatorial surface. For al l the models, checking the relativistic virial identities (jBonazzola fc GourgoulhonI Il994l : 
[Gci urgou lhon fc Bonazzolal Il994[ ). we confirm the typical values are orders of magnitude 10~^ — 10~^, which are 
acceptable values for present numerical scheme. 

3.2. Equations of State 

As mentioned in Sec. (TJ equation of state (EOS) is an important ingredient for determining the equilibrium configu- 
rations. Instead of conducting an extensive study as done before for the studies of rotating equi l ibrium configurations 
in which a greater variety of EOSs were employed, (e. g., iNozawa et aTl (1998); Morrisoii et al.l (120041) and references 
therein), we adopt here four k inds of EOSs, of SLy (|Do uchin fc Haensel 2001), FPS ( Pandharipande fc Ravenhalll 



Il98£i) . Shen (jShen et al.lll998ah . and Lattimer-Swesty(LS) (Lattimcr & Swcsty 1 99ll ). which are often employed in the 
recent MHD studies relevant for magnetars. 

In the study of cold neutron stars, the /3-equilibrium condition with respect to beta decays of the form + p < — > 
n + i/e and n < — > p+e~ +Pe, can be well validated as for the static properties. Since neutrinos and antineutrinos escape 
from the star, their chemical potentials vanish at zero-temperature T = with T being the temperature. Thus the /3- 
cquilibrium condition can be expressed as /in = Mc + /^p, with /io, and fip being the chemical potentials of neutron, 
electron, and proton, respectively. With the charge neutrality condition, we can determine the three independent 
thermodynamic variables, (for example the pressure as P(p,Ye,T), with being the el ectron fraction), can only b e 
determined by a single variable, which we take to be the density, namely P{p^Ye{p)) (jShapiro fc Teukolskvl fl983f l . 
noting here that T = is assumed for the case of the cold neutron stars. Thanks to this, we can use the formalisin 
mentioned in Sec. [5] without violating the barotropic condition of the EOSs. It should be noted that Rei seneggerl ()2008f ) 
has recently argued the breaking of the barotropic condition leads to the stabilization of the magnetic equilibria. This 
tells us the import ance to model the neutron stars by multicomponent fluid, such as neutrons, electrons, and protons 

eisenegge"ill2001[ ). However, taking into account this effect to equilibrium configuration is beyond scope of this article. 
Furthermore, the barotropic condition plays a crucial role to obtain the Bernoulli equation Eq. ([3|). Hence, in this 
article, we only consider the barotropic case for simplicity. 

At the maximum densities h igher than ~ 2pnuc, muons can appear, and higher than ^ '^Pm'- (jWiringa et al.lfl988l : 
lAkmal fc Pahdharipahdil998D . one may take into the possible appearance of hyperons (jGlendennin dl200lD. However, 



since the muon contribution to pressure at the higher density has been pointed to be verv small (jPouchin fc Haensell 
|2001| ). and we still do not have detailed knowledge of the hyperon interactions, we prefer to employ the above neutron 
star matter, namely e~,n,p, model to higher densities. In the following, we shortly summarize features of the EOSs 

employed here. 

The Lattimer-Swesty EOS (jLattimer fc Swestvlll991h has been used for many years as a standard in the research field 
of cor e- collapse supernova e xplosions and the subsequent neutron star formations (see references in ISumivoshi et all 
()2005() : iKotake et al.l (|2006f )): it is based on the compressible drop model for nuclei together with dripped nucleons. 
The values of nuclear parameters ar e chosen according to nuclear mass formulae and other theoretical studies with the 
Skyrme interaction. The Shen EOS (jShen et al.lll998ar ) is a rather modern one and currently often used in the research 
field. The Shen EOS is based on the relativistic mean field theory with a local density approximations, which has been 
constructed to rep roduce the experimental data of masses and radii of stable and unst able nuclei (see references in 
[^^FaDdillli)). FPS are modern version of an earlier microscopic EOS calculations bv lFriedman fc P andharipandd 
(Il981t) . which employs both two body (U14) and three-body interactions (TNI). In the SLy EOS (jPouc hin fc Haensel 
|2001|) . neutron-excess dependence is added to the FPS EOS, which i s more suitabl e for th e neutron star interiors. As 
for the FPS and SLy EOSs, we use the fitting formulae presented in lShibata et al.l ()2005[ ). 

4. NUMERICAL RESULTS 

In constructing one equilibrium sequence, we have three parameters to choose, namely the central density pc, the 
strength of the magnetic field parameter &, and the axis ratio rp/r^- Changing these parameters, we have to seek 
solutions in as wide parameter range as possible to study the properties of the equilibrium sequences in detail. For 
clarity, we categorize the comput ed models into two, namely, non-rotating or rotating and discuss them separately, as 
done in lKiuchi fc Yoshidal (|2008D . 

To calculate the non-rotating models, pc and 6 (see Eq. (|4])) have to be given with = 0, where Vl is the angular 
velocity. For the rotating models, one need specify the axis ratio rp/r^ in addition to pc and 6, where Tp and are 
the polar and the equatorial radii, respectively. To explore properties of the relativistic magnetized star models with 
the realistic EOSs, we construct thousands of the equilibrium models for each EOS, namely for the non-rotating cases, 
28 X 40 models in the parameter space of (pc, &), and for the rotating cases, 28 x 40 x 10 models in the parameter space 
of (pc h, rp/re). Throughout this section, we mainly pay attention to the magnetic field dominant models, i.e., the 
models with H > T, since we are concerned with the structure of the magnetar. For some magnetars, the condition of 
H > T could be realized because the internal magnetic field strength could reach an order of 10^* G and the typical 
spin period is several seconds. It should be noted that for the very fast spinning stars in which the dynamo process 
will work, the condition oi H < T might be realized and such models are important. However, recent observations 
support the fossil origin (jFerrario fc Wickramasinghell2006h rather than the dynamo process and our main interest is 
the observed magnetar so far. Hence, in this paper, we mainly refer to the models with H > T. 
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4.1. Non-rotating models 

First let us consider the static configurations for the following two reasons. (1) Since the magnetars and the high 
field neutron stars observed so far are all slow rotators, the static models could well be approximated to such stars. (2) 
In the static models, one can see purely magnetic effects on the equilibrium properties because there is no centrifugal 
force and all the stellar deformation is attributed to the magnetic stress. 

In Fig. [1] characteristic quantities are given for a M = 1.41M0 star with the SLy EOS characterized by pc = 
10^'^[g/cm^], Rcir = 14.4km, r^/re = 1.14, e = -0.75, = 0.187. In the figure, it can be seen that the strong 

toroidal magnetic fields peaking in the vicinity of the equatorial plane (panel (c)), acts through the Lorentz forces (see 
the arrows pointing perpendicular to the magnetic axis in panel (d)) to pinch the matter around the magnetic axis, 
making the stellar shape prolate (panel (a)). The maximum field strength reaches to 7.8 x lO^^G for this model. From 
panel (b), the regions where the magnetic pressure dominates over the matter pressure are seen to be formed. Indeed, 
the toroidal magnetic field lines behave like a rubber belt that is wrapped around the waist of the star. The matter 
fastened by this belt becomes stiff as seen from the large adiabatic indices shown in panel (e), where the definition 
of the adiabatic index is 91nP/91npQ|s with s being the entropy of the matter (s = for the cold neutron stars 
considered here). These gross properties are common to the other realistic EOSs. In the following, we move on to look 
the differences more in detail. 

Figures [2][4] show the equilibrium configurations for the other EOSs of FPS, Shen, and LS, plotted in the same style 
as Fig. [H The central density for each model is set to be the same. Comparing Fig. [T] to Fig. [21 the profiles of 
the stars with SLy and FPS EOSs are quite similar, reflecting the re semblance of the stiffness of the EOSs below 
Pc = 10^^[g/cm'^] considered here (See Fig. 1 in lKiuchi fc Kotcikg ()2008f )). On the other hand, the density distribution 
of the star with Shen EOS is found to become less prolate than those with SLy and FPS EOS (compare Fig. [TJa) with 
[3lja)). This can be understood by considering the differences in the distributions of the magnetic and matter pressures. 
The concentration of the magnetic field to the stellar center for Shen EOS is weaker than that for SLy or FPS EOS 
(compare Fig. [ijc) with[3jc)). Moreover, the stiffened matter with higher adiabatic indices extends further out for 
Shen EOS than for SLy or FPS EOS, (see Fig. [T]Je) and[3lje)), which means that the matter pressure stays relatively 
large up to the stellar surface for Shen EOS. As a result, the regions in which the magnetic pressure is dominant over 
the matter pressure can appear rather in the outer regions for Shen EOS than for SLy or FPS EOS (see Fig. [T]Jb) 
-IHlJc)). This implies that the magnetic fields for Shen EOS are effectively less fastening to pinch the matter around the 
magnetic axis than those for SLy or FPS EOS. For LS EOS, the regions in which the ratio of the magnetic pressure 
to the matter pressure, Pmag/^niattor, is large also exist near the stellar surface (see Fig. lUb)). However, the density 
distribution is found to become similar to that for SLy and FPS EOSs but not for Shen EOS because the magnitude 
of -Pmag/^matter IS Sufficiently higher than that for Shen EOS, producing the strong fastening (see the color legend of 
Fig. lib) and lib)). 

In Fig. [5l the gravitational mass A/ for the different EOSs with the constant magnetic flux sequences is shown as 
a function of the circumferential radius i?cir- Note that $30 means the flux normalized by units of lO'^^G cm^. and 
that the curves labeled by their values of $30 indicate the sequences with the constant magnetic flux. In each panel, 
the left thick red line represents the spherical star limit and the right thick red line does the non -convergence limit 
(b_2±lO in Eq. ([4])) beyond which any solution cannot converge with the present numerical scheme (jKiuchi fc Yoshidal 
l2008f l. Strictly speaking, values of b at the non-convergence limit depend on the central density pc- However, we 
confirm the non-convergence limit appears around 6 = 10 for all the central densities we set in this study. Thus, the 
non-convergence limits (the right thick red lines in Fig. [5]) look like smooth curves. For all the EOSs, the physically 
acceptable solutions exist in the regions bounded by these two lines. The filled circles on each line represents the 
maximum mass models. To focus on the sequence with the maximum mass, which we call as the maximum mass 
sequences henceforth, in Tabled we summarize their global physical quantities: the gravitational mass M, the baryon 
rest mass Mq, the circumferential radius at the equator i?cir, the maximum strength of the magnetic fields iJmax.is 
normalized by 10^* [G], the ratio of the magnetic energy to the gravitational energy lf/|VF|, and the mean deformation 
rate e. It should be noted that, irrespective of the EOSs, the sequences with $30 > 2.5 encounter the non-convergence 
limits before it reaches the maximum mass point (see also Table [1} . 

From Table [1] and Fig. O we can see that regardless of the different EOSs, the maximum gravitational masses 
M, do not increase with the magnetic fiux so much, in comparison with the other global quantities. This is in 
contr ast to the polytropic E OS, in which the maximum masses increase almost monotonically with the magnetic 
fiux (|Kiuchi fc Yoshidal l2008f ). For clarifying the relation, we prepare Fig. IHa), in which the increasing rates of the 
maximum masses, (M — AfsphcricaOZ-^'^sphcricab with A/gphcricai being the maximum masses of the spherical stars, are 
shown for the different EOSs along the maximum mass sequences. It is seen that the maximum masses, irrespective 
of the different EOSs, once decrease, and then begin to increase as i//|T4^| becomes larger. This first decrease is 
because the volume of the high density regions can temporary become smaller as the magnetized stars transform to 
the prolate shape due to the magnetic fastening, while keeping the changes in the central density small. In fact, the 
central densities of the maximum mass sequences stay almost constant with those of the spherical cases (see Table |T]) . 
For the polytropic EOS with n = 1 (referred as Pol2 in Fig. Hl^a)), with n being the polytropic constant, the different 
trend can be seen. The masses keep increasing as the function of i//|iy|. As shown in Table [1] (sec also Fig. |6|a)), 
the maximum gravitational mass of the neutron stars with $30 = 2.5 is larger than that with $30 = for FPS and 
LS EOS, but smaller for SLy and Shen EOS. It should be noted this does not imply the maximum mass with SLy or 
Shen EOS never increases because this is set by the non-convergence limit mentioned above. 
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Table 1 

Global physical quantities for the maximum gravitational mass models of the constant magnetic flux 

sequences of the non- rotating stars. 

^30 Pc.l5 M[Mq\ Mo[Mq\ flcir[km] g^ax.lS H/\W\ g 

SLy 

O.OOE+00 2.00E+00 2.05E+00 2.43E+00 9.96E+00 O.OOE+00 O.OOE+00 O.OOE+00 
1.50E+00 2.10E+00 2.02E+00 2.33E+00 1.03E+01 1.16E+00 8.88E-02 -1.71E-01 
2.00E+00 2.10E+00 2.01E+00 2.28E+00 1.09E+01 1.47E+00 1.47E-01 -3.14E-01 
2.50E+00 2.10E+00 2.02E+00 2.25E+00 1.17E+01 1.73E+00 2.10E-01 -5.03E-01 

EPS 

O.OOE+00 2.40E+00 1.80E+00 2.11E+00 9.27E+00 O.OOE+00 O.OOE+00 O.OOE+00 
1.50E+00 2.60E+00 1.78E+00 2.01E+00 9.89E+00 1.40E+00 1.17E-01 -2.48E-01 
2.00E+00 2.60E+00 1.78E+00 1.98E+00 1.07E+01 1.73E+00 1.86E-01 -4.53E-01 
2.50E+00 2.50E+00 1.81E+00 1.97E+00 1.22E+01 1.93E+00 2.61E-01 -7.55E-01 

Shell 

O.OOE+00 1.20E+00 2.42E+00 2.81E+00 1.34E+01 O.OOE+00 O.OOE+00 O.OOE+00 
1.50E+00 1.30E+00 2.39E+00 2.72E+00 1.39E+01 7.26E-01 7.83E-02 -1.71E-01 
2.00E+00 1.30E+00 2.39E+00 2.69E+00 1.46E+01 9.09E-01 1.29E-01 -3.11E-01 
2.50E+00 1.30E+00 2.40E+00 2.66E+00 1.57E+01 1.06E+00 1.82E-01 -4.95E-01 

LS 

O.OOE+00 2.00E+00 1.93E+00 2.23E+00 1.04E+01 O.OOE+00 O.OOE+00 O.OOE+00 
1.50E+00 2.10E+00 1.92E+00 2.15E+00 1.12E+01 1.16E+00 1.07E-01 -2.31E-01 
2.00E+00 2.10E+00 1.93E+00 2.13E+00 1.21E+01 1.43E+00 1.70E-01 -4.18E-01 
2.50E+00 2.00E+00 1.96E+00 2.13E+00 1.37E+01 1.58E+00 2.39E-01 -6.98E-01 



Fig. [DJb) shows the degree of the deformation for the different EOSs as a function of before the non- 

convergence Hmit. It is found that the degree can reach to an order of e = — O(O.l) irrespective of the EOSs. 
Comparing with the polytropic case (Pol2), we find the mean deformation rates become smaller for the models with 
the realistic EOSs. The reason is as follows. The magnetic belt effects, due to which the matter is pinched, subside in 
the vicinity near the equatorial plane less than ^ 10 km (sec Fig.[TId) -Hl^d)). The adiabatic indices there are generally 
higher than 2 (see Fig. [He)-|He)) irrespective of the different EOSs. This means that the matter is stiffer than the 
case of n = 1 polytropic stars. As a result, the magnetic fastening becomes less, leading to the smaller deformation 
for the realistic EOSs. 

4.2. Rotating models 

Next, we proceed to the models with rotation. In Fig. [71 we display typical distribution of the rotating equilibrium 
star taking a M = 2.00Mo model with the SLy EOS characterized by pc ~ 6.0 x 10^^[g/cm'^], i?cir = 25.0km, 
rp/re ^ 0.68, e = -0.24, T/\W\ = 4.77 x 10"^ Hl\W\ = 0.152. This model has the maximum magnetic field strength 
of 5.49 X lO^^G and rotates in the mass-shedding limit. Comparing Fig. [7| with Fig. [1] we see that basic properties of 
the rotational equilibrium structure are similar to those of the static models. Difference between them only appears 
near the equatorial surface of the stars. As shown in Fig. [71 the density distribution is stretched from the rotation axis 
outward due to the centrifugal force and the shape of the stellar surface becomes oblate, which can be also seen from 
the value of the axis ratio Tp/ve- On the other hand, the equi-density contours deep inside the stars are prolate, as 
confirmed from the value of the mean deformation rate e. In the central regions, the magnetic fastening of the matter 
generically dominates over the centrifugal forces. 

It is should be emphasized that we here single out the magnetic field dominated models in which the ratio of the 
rotational energy to the gravitational energy T/|VF| is much smaller than the ratio of the magnetic energy to the 
gravitational energy even when they rotate at the breakup angular velocity because we focus on the magnetic 

effects on the stellar structure. In fact, typical values oi H/\W\ and T/|VF| amongst the models computed in this study 
are orders of 10^^ and 10^^, respectively. (As we will see later, we cannot pick out such magnetic field dominated 
models for the supramassive sequences of the magnetized stars with the realistic EOSs.) The effects of the rotation 
on the stellar structures are therefore secondary in the models shown in this subsection. Likewise, the effects of the 
EOS on the rotational equilibrium sequences are predominantly determined by the ones in the non-rotating sequences, 
which we mentioned in the previous section. 

4.3. Constant baryon mass sequence 

In this su bsection, we pay a t tentio n to the constant baryon mass sequences of the rotating magnetized stars. As 
discussed in iKiuchi fc Yoshidal ()2008f ). it is possible to single out a sequence of the equilibrium stars by keeping the 
baryon rest mass and the magnetic flux constant simultaneously. Such equilibrium sequences may model the isolated 
neutron stars that evolve adiabatically losing angular momentum via the gravitational/electromagnetic radiations. For 
simplicity, we omit the evolution in the function K (Eq. ([2])), which will change as a result of losing angular momentum, 
although there has been no simple way to determine it so far. It should be not ed that the efficient gravitational 
radiations can be made possible only if the spin and magnetic axes are misaligned (ICutleij|2002l) . However, including 
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the misalignment to our model is extremely difficult because it breaks down the equilibrium condition due to the 
gravitational radiation. With the al ignment axes, the spin-do wn via electromagnetic radiations could be possible if the 
outside of the star is non- vacuum (jGoldreich fc Juliaiilll969[ ). However, the condition that the outside of the star is 
vacuum is crucial here for calculating the equilibrium configurations. Furthermore, the method to take into account the 
mixed (poloidal-toroidal) field in general relativistic equilibrium configuration is not still established as we mentioned 
in Sec. [T] Bearing these caveats, we boldly pay attention to the constant baryon mass and magnetic flux sequence 
obtained fro m the equilibrium configurations to mimic the evolution of the magnetized neutron star in the following. 

As done in lCook et al.l (jl992[ ). we divide the equilibrium sequences of the magnetized stars into two classes, normal 
and supramassive equilibrium sequences. In this study, the normal (supramassive) equilibrium sequence is defined as 
an equilibrium sequence whose baryon rest mass is smaller (larger) than the maximum baryon rest mass of the non- 
magnetized and non-rotating stars. In Tabled the maximum baryon rest mass of the spherical (magnetized rotating) 
stars, which are referred as Mo,sph,max(Afo,max), with the each EOSs, are given. Note that Afo,max is the maximum 
baryon mass of the magnetized star models calculated in this study. 

Table 2 

Maximum baryon mass of the spherical stars and the magnetized rotating star with the each EOSs. 



SLy FPS Shen LS 

Mo,sph,max[M0] 2.43 2.11 2.81 2.23 
Mo,inax[MQ] 2.82 2.46 3.36 2.62 

4.3.1. Normal sequence 

First, let us consider the normal equilibrium sequences of the magnetized rotating stars, choosing Mq ~ 1.90Mq 
stars with the magnetic flux of $30 = 1. Here Table [2] tells us that all the stars with Mo = 1.90Mq belong to the 
normal sequence irrespective of the EOS. The global physical quantities for the models with the different EOSs are 
summarized in Table [3] for convenience. It should be emphasized that we select the magnetic field-dominant sequence 
(ff/|W^| > r/|M^|) because we are interested in the effect of strong magnetic field on the stellar structure (see Table [3]). 

In Fig.m the relative change of the global physical quantities to those of the spherical stars, such as M, i?cir, Pc,i5, 
J, ^max.iSi a-iid e are given as functions of for the constant baryon mass and magnetic flux equilibrium sequences 
for the different EOSs. Here fl^ means the angular velocity normalized by 10^[rad/s]. It can be seen from panel 
(a), all the normal equilibrium sequences, starting at the non-rotating equilibrium stars = 0), continue to the 
mass-shedding limits, at which the sequences terminate (note in each panel, the asterisks indicate the mass-shedding 
limits). As the angular velocity increases, the gravitational masses and the circumferential radii increase (panels (a) 
and (b)) due to the centrifugal forces, while the central densities and the maximum amplitudes of the magnetic fields 
decrease (panels (c) and (e)) because of the constancy of the baryon rest mass and the magnetic flux, respectively. 
The angular momenta monotonically increase with the angular velocity, which implies that the angular momentum 
loss via gravitational radiation results in the spin down of the star for the normal equilibrium sequences. As seen, 
these qualitative behaviors are found to be common to the realistic EOSs employed in this study. 

Choosing SLy EOS for example, we show in Fig.[5]how the masses and the angular momenta change with the angular 
velocity along with the constant magnetic flux. In each panel, the non-magnetized case ($30 = 0) is given for the 
sake of comparison. As seen, qualitative changes induced by rotation are basically independent of the value of $39. 
We find that the angular velocity at the mass-shedding point becomes small due to the Lorenz force as the magnetic 
flux increases. Figure[7{d) (or Figure[T](d)) clearly shows that near the stellar surface, the Lorenz force exerted on the 
matter has the same direction with the centrifugal force. Therefore, the stars with the strong magnetic flelds are easy 
to reach the mass-shedding limit. 

To see clearly the dependence of the mass-shedding angular velocity fims on the baryon mass and magnetic flux 
characterizing the equilibrium sequences, we draw the equi-fims contours plot on Mo-$ plane in Fig. 1101 It can be 
seen that irrespective of the EOSs, a large magnetic flux decreases substantially the mass-shedding angular velocities. 
Comparing ri„is for the normal equilibrium sequences with the same baryon mass and magnetic flux, we also find that 
for the Shen EOS models Jims takes the smallest values, for the LS and SLy EOS models middle values, and for the 
FPS EOS models the largest values (see also Fig.[8ja)). This implies we may constrain the EOSs via angular velocity 
observations. For example, if we observe the magnetized neutron star with Mq ~ 2.0Mq and = 8 x lO'^rad/s, Shen 
EOS is rejected because the stars with Shen EOS cannot rotate at 17 = 8 x lO'^rad/s irrespective of the strength of the 
interior magnetic fields. 

4.3.2. Supramassive sequence 

Now we move on to consider the supramassive equilibrium sequences of the magnetized rotating stars. As in the 
case of the normal sequence, the qualitative features in the sequences are not sensitive to the difference in the EOSs. 
So we take the model with Shen EOS in the following. Figure [TT] gives the plots of the global stellar quantities as a 
function of ^l for a constant mass (isupramassive Mq = 2.9OM0) and magnetic flux with the Shen EOS. Each curve is 
labeled by its value of <I>3o which is held constant along the sequence. Values of the global physical quantities, pc, M, 
i?cir, ^, J, T/\W\, H/\W\, e, and Bmax.iS: for some selected models are summarized for convenience in TablelH 
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Table 3 

Global physical quantities for the normal equilibrium sequences of the rotating stars with Mq — 1.90Mq 

AND $30 = 1. 



Pc,15 M[Mq\ flcirlkm] 


^3 


cJ/GM'^ T/\W\ H/\W\ 


e 


-Bmax.lS 


SLy 



1.21E+00 1.70E+00 1.21E+01 O.OOE+00 O.OOE+00 O.OOE+00 9.40E-02 -2.59E-01 7.33E-01 

1.17E+00 1.71E+00 1.28E+01 4.25E+00 7.91E-01 2.27E-02 9.80E-02 -1.70E-01 7.29E-01 

1.13E+00 1.72E+00 1.44E+01 5.97E+00 1.21E+00 5.01E-02 1.03E-01 -7.09E-02 7.19E-01 

l.llE+00 1.72E+00 1.57E+01 6.38E+00 1.33E+00 6.00E-02 1.04E-01 -3.61E-02 7.14E-01 

l.lOE+00 1.73E+00 1.73E+01 6.50E+00 1.37E+00 6.33E-02 1.05E-01 -2.40E-02 7.12E-01 



MS - - - - - - - - 

FPS 

1.75E+00 1.68E+00 1.07E+01 O.OOE+00 O.OOE+00 O.OOE+00 8.51E-02 -2.09E-01 8.82E-01 
1.61E+00 1.70E+00 1.18E+01 5.88E+00 9.54E-01 3.38E-02 9.27E-02 -9.65E-02 8.61E-01 
1.57E+00 1.70E+00 1.23E+01 6.53E+00 1.09E+00 4.36E-02 9.48E-02 -6.55E-02 8.54E-01 
1.52E+00 1.71E+00 1.35E+01 7.30E+00 1.29E+00 5.85E-02 9.80E-02 -1.89E-02 8.41E-01 
1.49E+00 1.71E+00 1.56E+01 7.62E+00 1.38E+00 6.67E-02 9.97E-02 7.00E-03 8.34E-01 

MS - - - - - - - - 

Shen 

5.47E-01 1.75E+00 1.66E+01 O.OOE+00 O.OOE+00 O.OOE+00 1.04E-01 -3.55E-01 4.22E-01 
5.32E-01 1.75E+00 1.77E+01 2.61E+00 8.18E-01 2.03E-02 1.07E-01 -2.58E-01 4.19E-01 
5.18E-01 1.76E+00 1.91E+01 3.44E+00 1.14E+00 3.76E-02 l.lOE-01 -1.80E-01 4.15E-01 
5.08E-01 1.76E+00 2.08E+01 3.85E+00 1.32E+00 4.97E-02 l.llE-01 -1.28E-01 4.12E-01 
5.03E-01 1.77E+00 2.29E+01 4.02E+00 1.41E+00 5.59E-02 1.12E-01 -1.02E-01 4.10E-01 

MS - - - - - - - - 

LS 

1.18E+00 1.71E+00 1.28E+01 O.OOE+00 O.OOE+00 O.OOE+00 9.14E-02 -2.55E-01 6.65E-01 
1.15E+00 1.72E+00 1.33E+01 2.79E+00 5.61E-01 l.llE-02 9.35E-02 -2.10E-01 6.59E-01 
1.07E+00 1.73E+00 1.49E+01 4.99E+00 l.llE+00 4.12E-02 9.92E-02 -9.80E-02 6.37E-01 
1.03E+00 1.73E+00 1.63E+01 5.54E+00 1.30E+00 5.47E-02 1.02E-01 -5.02E-02 6.26E-01 
l.OlE+00 1.74E+00 1.88E+01 5.78E+00 1.40E+00 6.23E-02 1.03E-01 -2.37E-02 6.19E-01 

MS - - - - _ _ _ _ 



From panel (a), the supramassive sequences begin at the mass-shedding hmits with lower central density pc as the 
angular velocity decreases and reach the turning points (the slowest rotation points), then move to the other mass- 
shedding limits with higher central density pc as the angular velocity increases. The gravitational mass decreases at 
first, and encounters the turning point, then it begins to increase as the central densities decrease (panel (b)). The 
circumferential radius decreases because of the increase of the central density and keeps the almost constant value 
after the turning point (panel (c)). From Figs. [TlT d). we find most models of the supramassive sequences have the 
positive values of the mean deformation rate, which implies the stars oblately deform. This feature can be also seen 
in the sequences with the other EOSs. 

In fact, Table |4] shows that the rotational energy dominates over the ma gnetic energy (T/lW'l > in these 

seque nces. This means that we pick out the rotation-dominated sequence ( Yoshida fc Eriguchil l2006l : lYoshida et al.l 
[200l for Table 4 and Fi gure 11. However, it should be noted that it depends on the values of the baryon rest mass 
and magnetic flux whether a sequence belongs to the rotation-dominated one or the magnetic-field-dominated one. To 
get an overall picture, we show the phase diagram in Figure [T^] on Mq and $30 plane of the supramassive sequence 
with eac h realistic EOS. Th e solid red line in Fig. [T^] indicates the boundary, above which there is no physical solution 
(see, e.g. lKiuchi fc Yoshidal (|2008[ )). The regions below the solid lines correspond to a supramassive sequence. We 
plot the sequences given in Table U as the cross symbols in each panel. In the right side regions of the dotted lines in 
Fig. [T2I the sequences only have the oblate (e > 0) stars. On the other hand, in the left side regions, there exit the 
sequences which have the prolate (e < 0) stars (see also the magnified panels in the each plots). Figure [T^] clearly shows 
that the supramassive sequence with the realistic EOSs almost belong to the rotation-dominated sequence because 
the magnetic-field-dominated sequences are easy to be prolately deformed as discussed in Sec. 14.21 Looking carefully, 
it can be seen that for FPS EOS, the parameter region which permits the prolately deformed stars is widest. For 
Lattimer-Swesty EOS, the region is narrowest vice versa and for SLy and Shen EOS, it is in medium. It should be 
mentioned th at for the poly trop i c EO S, the magnetic-field-dominated sequences along the supramassive sequences 
easily appear (jKiuchi fc Yoshidal (|2008f l) and this difference may be useful for the constraining EOS via gravitational 
wave observations as we will discuss in Sec. [SJ 

Finally, let us investigate the spin up of the stars as the stellar angular momentum decreases, which is one of the 
astrophysically important features o f the co nstant baryon m ass equilibrium sequences of the non-magnetized rotating 
stars (see, e.g. IShapiro et all (|1990f ) and iCook et al.l ()1992f )). This spin up effect of the rela tivistic stars con t aining 
pure ly poloidal magne t ic fiel ds and purely toroidal magnetic fields has also been found by iBoucauet et al.l (|1995f ) 
and iKiuchi fc Yoshidal ()2008D . respectively. In this study, we find that the the spin up effects occur for the stars 
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Table 4 

Global physical quantities for the supramassive equilibrium sequences of the 

rotating stars. 



Pc, 15 


M[Mq] iJcir[km] ^3 cJ/GA/| TI\W\ UI\W\ 


e 




SLy (A/o,<I>3o) = {2.45Af0,l.O) 



MS 



1.27E+00 2.14E+00 1.51E+01 8.95E+00 2.89E+00 l.lOE-01 2.02E-02 2.26E-01 3.93E-01 
1.50E+00 2.11E+00 1.18E+01 7.69E+00 2.08E+00 6.14E-02 1.63E-02 1.26E-01 4.01E-01 
2.00E+00 2.08E+00 1.02E+01 4.67E+00 1.04E+00 1.65E-02 1.15E-02 2.39E-02 4.11E-01 
2.60E+00 2.11E+00 9.79E+00 9.29E+00 1.96E+00 4.95E-02 9.28E-03 1.04E-01 4.30E-01 
2.80E+00 2.13E+00 9.85E+00 l.llE+01 2.35E+00 6.83E-02 8.68E-03 1.46E-01 4.35E-01 
3.00E+00 2.15E+00 l.OOE+01 1.27E+01 2.70E+00 8.64E-02 8.10E-03 1.86E-01 4.39E-01 

MS - - - - - - - - 

EPS (Afo,^3o) = (2.2OAf0,l.O) 

MS - - - _____ 

1.58E+00 1.93E+00 1.40E+01 9.52E+00 2.29E+00 1.04E-01 2.38E-02 2.13E-01 4.60E-01 
1.80E+00 1.92E+00 1.15E+01 8.97E+00 1.91E+00 7.53E-02 2.08E-02 1.55E-01 4.74E-01 
2.30E+00 1.90E+00 l.OlE+01 8.20E+00 1.48E+00 4.65E-02 1.62E-02 9.31E-02 4.99E-01 
2.60E+00 1.90E+00 9.74E+00 8.84E+00 1.52E+00 4.81E-02 1.46E-02 9.90E-02 5.13E-01 
2.80E+00 1.90E+00 9.62E+00 9.69E+00 1.64E+00 5.44E-02 1.37E-02 1.15E-01 5.22E-01 
3.00E+00 1.91E+00 9.57E+00 1.08E+01 1.79E+00 6.34E-02 1.30E-02 1.36E-01 5.31E-01 

MS - - - - - - - - 

Shen (Mo,'I>3o) = (2.9AfQ,1.0) 

MS - - - _____ 

7.78E-01 2.57E+00 2.05E+01 6.21E+00 3.80E+00 9.32E-02 5.16E-02 1.59E-01 4.38E-01 
8.04E-01 2.56E+00 1.82E+01 6.08E+00 3.60E+00 8.47E-02 5.00E-02 1.40E-01 4.43E-01 
1.19E+00 2.53E+00 1.46E+01 5.16E+00 2.40E+00 4.02E-02 3.78E-02 4.59E-02 4.92E-01 
1.40E+00 2.53E+00 1.42E+01 5.80E+00 2.58E+00 4.51E-02 3.47E-02 6.72E-02 5.14E-01 
1.80E+00 2.56E+00 1.44E+01 7.80E+00 3.38E+00 7.18E-02 3.11E-02 1.46E-01 5.53E-01 
2.00E+00 2.58E+00 1.52E+01 8.70E+00 3.78E+00 8.65E-02 2.99E-02 1.88E-01 5.69E-01 

MS - - - - - - - - 

LS (A/q.'I'so) = (2.3A/0,1.O) 

MS - - - _____ 

1.14E+00 2.05E+00 1.66E+01 7.65E+00 2.58E+00 1.02E-01 2.29E-02 2.19E-01 3.52E-01 
1.30E+00 2.03E+00 1.35E+01 7.24E+00 2.13E+00 7.38E-02 2.03E-02 1.59E-01 3.70E-01 
1.80E+00 2.01E+00 1.13E+01 5.95E+00 1.38E+00 3.27E-02 1.52E-02 6.42E-02 4.10E-01 
2.00E+00 2.00E+00 l.lOE+01 6.14E+00 1.35E+00 3.11E-02 1.40E-02 6.22E-02 4.24E-01 
2.30E+00 2.01E+00 1.07E+01 7.32E+00 1.54E+00 3.92E-02 1.26E-02 8.55E-02 4.44E-01 
2.40E+00 2.01E+00 1.06E+01 7.91E+00 1.65E+00 4.41E-02 1.22E-02 9.85E-02 4.50E-01 

MS - - - _____ 



containing purely toroidal magnetic fields with the realistic equations of state. Figure fTSl a) shows the angular velocity 
17 as a function of the angular momentum J along the supramassive sequences with LS EOS. These sequences are 
characterized by the constant magnetic flux $30 = 1.0 and the number attached by each curve indicates the constant 
baryon mass values in units of Mq. We observe the spin up effect, the angular velocity increase as the angular 
momentum decreases, in the supramassive sequences with 2.39, 2.41, and 2.43AfQ. This spin up is induced by the 
decrease of the moment of inertia with decreasing the angular momentum. Reaching the turning point of the constant 
sequ ence path of Fig. 13, the stars in absent of th e magnetic fields were pointed out to begin to collapse to black 
hole (jCook et al.l (|1992t ) according to the criterion bv lFriedman et aTl (|1988l )). Qualitatively the magnetized stars here 
are also expected to form black holes on the analogy, but to find the criterion for the magnetized case needs further 
investigation. We here define a critical angular velocities, f2up, with which the stars can start to spin up. Finally, 
we get the relationships between fiup and $ for the different EOSs, which are drawn in Fig. [TSTb). Our analysis 
indicates that the hierarchy of this spin-up angular velocity is l^up.SLy > f^up,FPS > ^^up.LS > ^^up.shen and this relation 
holds even if the sequences have strong magnetic fields. Furthermore the values of fiup are found to increase with the 
constant magnetic flux irrespective of the EOSs and the lines for the different EOSs never cross each other. 

5. DISCUSSION AND SUMMARY 
5.1. Discussion 

We give discussions in this subsection, paying attention to the properties of the strongly magnetized neutrons stars 
with the realistic EOSs. 

(1) As mentioned, the magnetized neutron stars could evolve adiab atically losing its angular momentum by grav- 
itational or electromagnetic radiation. Equations (2.9) and (2.10) in ICutled (|2002| ) give us the order estimation of 
angular momentum loss time scale by gravitational radiation tqw and that by electromagnetic radiation tem as 
'''Gw/'''EM = 8.9 X 10^^(i3d/10^G)^(eB/10^'')^^(i^s/kHz)^^, where e^, Bd, and I's represent the degree of deformation 
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induced by the toroidal magnetic field, the dipole magnetic field outside the star, and the rotation frequency. If we 
assume Bd = lO^^G and es = 10^^, the equation implies the angular momentum loss is driven by the electromagnetic 
radiation. On the other hand, our models suggest that eb can reach to the order of 0.1 as we discussed in Sec. 4. In 
this case, the gravitational radiation could be main agent of the angular momentum loss. 

As mentioned, the mean deformation rates e for the strongly magnetized stars depend on the parameters of the 
equilibrium models. Along the normal sequences, they are basically negative, which means that the mean matter 
distributions are prolate, even when the stars rotate at nearly the breakup angular velocity (see panel (f) of Fig. [5]). 
High toroidal fields enough to make the equilibri um configura tion prolate in the absence of rotation, are pointed out 
to be a good emitter of the gravitational waves (|Cutlerj|2002l ). leading to a secular instability, in which the wobbling 
angles between the rotation axis and the star's magnetic axis would grow on the dissipation timescale, until they 
become orthogonal. 

Following to Eq. (4.2) in loltfal (j2002), S/N = 11.7 (lOkpc/i^) (es/lQ-e) {W^iQ/Bd) [In (/,„ax//min)], we may 
estimate the signal to noise ratio of gravitational wave from the magnetized neutron star, where D and /max(/min) 
represent the distance to the source and the maximum (minimum) rotation frequency. Here, it is assumed that the 
neutron star spin-downs due to the magnetic dipole radiation, though in the present neutron star models, we fully 
omit the exterior dipole magnetic field B^ assuming it does not affect the neutron star structures. In our models, we 
find the maximum value of is of an order of 10~^, if the magnetic fields nearly reach the equipartition. Supposing 
the dipole magnetic field is lO^^'G, we have S/N 10 for D = 1 Gpc. To constrain the EOSs by the gravitational wave 
observations, we need to investigate the gravitational waveforms in details using the presented models here, which will 
be presented soon elsewhere. 

(2) The stars in the supramassive sequences with the sufficient large baryon mass start to spin up losing their 
angular momentum, irrespective of the EOSs. The critical angular velocities at which the spin-up begins depend on 
the specified constant magnetic flux and the selected EOS. Let us discuss the possibility of constraining the equation 
of state by the spin-up effect. If we detect the spin-up effect in a magnetized neutron star and estimate its critical 
angular velocity as e.g., 8 x 10'^[rad/s], SLy and FPS EOS would be rejected because the spin- up never occurs in the 
stars with these two EOSs. In our models which can spin up losing their angular momentum, the ratio of the rotational 
energy to the gravitational biding energy T/|I^| never excess 0.27, which is the conservative value of the dynamical 
instability. This result implies the strongly magnetized stars may survive during its spinning up. 

(3) The gravitational masses of the strongly magnetized stars constructed in this study are seemingly too large (see 
Ta ble IH because the canonical value of the neutron star mass is 1.4M0. However, the population synthesis study 
in iFerrario fc Wickramasinghd ()2006f ) indicates the neutron stars can possess the heavy mass ~ 2Mq if they have 
the strong magnetic fields. Thus, our strongly magnetized neutron stars could be models of the magnetar and the 
high field neutron star. In combination with (1), we may obtain the information about the field strength from the 
observation of the gravitational waves. 

(4) We shall comment on the definition of the su pramassive a nd th e normal sequences for the magnetized stars, 
which are categorized by the baryon mass following ICook et al.l ()1992f ) in this study. For the non-magnetized stars, 
a characteristic feature of the normal sequences is that they always begin at a non-rotating solution and end at a 
mass-shedding solution. For the magnetized stars, on the other hand, we find that some normal sequences do not 
start at a non-rotating solution and include no non-rotating solutions. This implies that for the magnetized stars, 
we have yet an another option to define the supramassive and the normal sequences. The definition can be given by 
the way how the equilibrium sequences end. If we define the supramassive sequences as the sequence that terminates 
at the mass-shedding solutions, we could extend the supramassive sequences to the lower mass regions (see Fig. I12p . 
As can be expected from Fig. [121 these extended supramassive sequences allow more prolate solutions, which should 
be inevitably magnetic-field dominated. However seeking the solutions in such parameter regime is hindered due to 
the non-convergence limits. We are now undertaking to update our numerical scheme to make the phase diagram of 
Fig. [12] more c omprehensive . 

(5) Recently, iKiuchi et al.l (|2008f ) have performed axisymmetric stability analysis of the toroidal magnetic field by 
making the GRMHD simulation and found the magnetic field with 7^ 1 is unstable. Therefore, we limit ourselves to 
fc = 1 case in this work. As for non-axisymmetric stability, Tayler has shown the toroidal magnetic field with fc = 1 
induces the kink instability near the magnetic axis in the Newtonian analysis l|Tavleilll973f ). Therefore, some models 
obtained in this study might be unstable due to the kink instability. However, the perturbation analysis only predicts 
the occurrence of instability. It is still unclear how the magnetic field evolves a fter the Tayler inst ability sets in. Will 
they decay or settle down to "new" equilibrium state like the axisymmetric case (iKiuchi et al.ll200"8l) ? To draw a robust 
conclusion, the 3D GR MHD simulations are necessary, which is our future work. 

(6) As mentioned in section 1.1, the strong magnetic fields inside neutron stars shoul d decay during t heir c volu- 
tions, which has been totally neglected in this study. Employing Eq. (61) in .Goldreic h fc Reiseneggerl ()1992| ). the 
lifetime (decaying timescales of the magnetic fields) for the presented models here can be estimated as an order of 
~ 10'^ years. Although more rapid decay is thought to be possible by taking into account the neutron star crusts 
([Pons fc Geppertll2007f). the lifetime s eems not so contradictory with the observations implying the lifetimes less than 
10'^ years ([Woods fc Thompsonll2004l ). 

(7) In this paper, we considered only the magnetic effects on the equilibrium configurations through the Lorentz force 
but did not take into account the additional changes caused by the magnetic effects on the EOS. For the super-strong 
magnetic fields of i? > lO^iJqED, there appear two counter effects in the EOSs, namely the stiffening and softening due 
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to th e anomalous magnetic moments of the nucleons and to the Landau quantization, respectively (jBroderick et al] 
Since the equilibrium configurations are determined by the subtle local balance of the pressure gradient, vs. 
predominantly the gravitational force (plus the Lorentz force and the centrifugal force), it is by no means a trivial 
problem to see the effects of the local change of the pressure on the important global properties of the strongly 
magnetized equilibrium star, such as the mass and the radius. To answer this important problem, it is indispensable 
to incorporate the magnetic corrections to the EOSs, however beyond scope of this paper. 

5.2. Summary 

In this study, we have investigated equilibrium sequences of relativistic stars containing purely toroidal magnetic fields 
with four kinds of realistic EOSs of SLy, FPS, Shen, and LS, which have been often employed in recent MHD studies 
relevant for magnetized neutron stars. Solving master equations numerically using the Kiuchi-Yoshida scheme, we have 
constructed thousands of equilibrium configurations in order to study the effects of the realistic EOSs. Particularly 
we have paid attention to the equilibrium sequences of constant baryon mass and/or constant magnetic flux, which 
model evolutions of an isolated neutron star, losing angular momentum via gravitational waves. Important properties 
obtained in this study are summarized as follows ; (1) Along the maximum gravitational mass sequences, the stars 
with the realistic EOSs cannot be prolate as much as the stars with a polytropic EOS of n = 1 with n being the 
polytropic index, which has been often employed to mimic the nuclear EOSs. (2) The dependence of the mass- 
shedding angular velocity on the EOSs along a constant baryon mass and magnetic flux sequences, is determined 
from that of the non-magnetized case. Along the sequence, the stars with Shen(FPS) EOS reach the mass-shedding 
limit at the smallest(largest) angular velocity, while the stars with SLy or Lattimer-Swesty EOSs take the moderate 
values. (3) For the supramassive sequences, the equilibrium configurations become generally oblate for the realistic 
EOSs, although the prolately deformed stars can exist in a narrow parameter region spanned by the constant baryon 
mass and magnetic flux. For FPS EOS, the parameter region which permits the prolately deformed stars is widest. 
For Lattimer-Swesty EOS, the region is narrowest vice versa and for SLy and Shen EOS, it is in medium. (4) For 
the supramassive sequences, the angular velocities fiup, above which the stars start to spin up as they lose angular 
momentum, are found to depend sharply on the realistic EOSs. The hierarchy of the critical spin up angular velocity 
is SLy > FPS > LS > Shen EOS and this turn never change even if they have strong magnetic fields. Our results 
suggest the relativistic stars containing purely toroidal magnetic fields will be a potential source of gravitational waves 
and the EOSs within such stars can be constrained by observing the angular velocities, the gravitational waves, and 
the signature of the spin up. 

K. Kiuchi expresses thanks to K. Maeda and S. Yamada for invaluable discussion. K.Kotake expresses thanks to 
K. Sato for continuing encouragements. Numerical computations were in part carried on XT4 and general common 
use computer system at the center for Computational Astrophysics, CfCA, the National Astronomical Observatory 
of Japan. This study was supported in part by the Grants-in-Aid for the Scientific Research from the Ministry of 
Education, Science and Culture of Japan (Nos. 19540309 and 20740150). 
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Fig. 5. — Circumferential radius versus gravitational mass for models with (a) SLy, (b) FPS, (c) Slien, and (d) LS EOS. The left thick 
line represents spherical, non-magnetized configurations, and the right thick line represents the boundary beyond which converged solution 
could not be obtained. Along the dotted lines labeled with >I>3o, the configurations have constant magnetic flux, and <i>3o denotes the flux 
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given in Table [T] 




0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 

H/|W| H/|W| 
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Fig. 8. — Global physical quantities Af, iJcin Pc, J, Smax, and e for the constant baryon mass and magnetic flux equilibrium sequences 
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SLy, FPS, Shen, LS EOS, respectively. All the plots show the increasing (decreasing) rates of the global quantities from those for the 
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Fig. 9. — Global physical quantities M and J for the constant baryon mass and magnetic flux equilibrium sequences with SLy EOS, 
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Fig. 12. — Phase diagram on constant baryon mass and magnetic flux plane of the supramassive sequences for (a)SLy, (b) FPS, (c) Shen, 
and (d) LS EOS. The solid lines correspond to the critical line, above which there exists no physical solution (see the text for details). The 
dashed lines show the models with e = 0. In the left (right) side regions of the dotted lines, the supramassive sequences have the prolate 
(oblate) solutions. The magnified panels arc put on the upper right of the each figures. The cross symbols in each panel represent the 
sequences given in Table |4] (see the text for details) . 




Fig. 13. — (a) Angular momentum and angular velocity relation in the supramassive sequences with LS EOS and <I>3o = 1.0. The 
numbers attached on the each sequences represent the values of constant baryon rest mass in units of Mq. (b) Critical spin up angular 
velocity as functions of the magnetic flux for each EOS. The critical angular velocity is defined by the value from which the stars start to 
spin up (see the text in details). 



